Thermopower of a superconducting single-electron transistor 
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Abstract 

We present a linear-response theory for the thermopower of a single-electron transistor consisting 
of a superconducting island weakly coupled to two normal-conducting leads (NSN SET). The 
thermopower shows oscillations with the same periodicity as the conductance and is rather sensitive 
to the size of the superconducting gap A. In particular, the previously studied sawtooth-like shape 
of the thermopower for a normal-conducting single-electron device is qualitatively changed even 
for small gap energies. 
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The transport properties of small conducting grains in the Coulomb blockade (CB) regime 
have extensively been studied during the past years. This regime is characterized by a new 
energy scale, the so-called charging energy E c of the grain (see below). The most prominent 
phenomenon is the occurence of CB oscillations in the low-temperature conductance of a 
small grain weakly coupled to the leads [l[ . Recently, thermoelectric effects in single-electron 

devices such as the thermopower have attracted growing interest DDDDDDQDD mi. 

The thermopower is related to the current that arises due to a finite temperature difference 
between the two leads It yields additional information about the kinetics of the system 
as it measures the average energy of the electrons carrying the current through the system. 
Therefore some type of electron-hole asymmetry in the system is necessary in order to 
observe a non-vanishing thermopower. 

In analogy to the CB oscillations of the conductance the thermopower of a small grain 
shows oscillations of the same periodicity but with sawtooth-like shape 

a a. 

In con- 
trast to the conductance this dependence on the external gate voltage is very sensitive to 
the conditions under which the thermoelectric transport occurs. This sensitivity has been 
demonstrated, e.g., for the transition from the sequential tunneling regime to the cotunnel- 
mg regime Recently, the thermopower of open quantum dots with strong coupling 

to the leads was investigated 0, Q, Further, the influence of Kondo correlations in 



ultra-small quantum dots on the thermoelectric effects was studied in Refs. |9|,llQ] while the 
thermopower of a molecule with internal degrees of freedom and weakly coupled to the leads 
was discussed in Ref. jlll ]. 

It is surprizing that, despite the enormous interest in superconducting SETs, the ther- 
mopower of such structures has not been investigated yet. In this work we study theoretically 
the thermopower of an NSN SET, i.e., a small superconducting island that is weakly coupled 
to normal-conducting leads (cf. Fig. [IJ, in an experimentally accessible regime. We show 
that even for rather small superconducting gaps (compared to the charging energy of the 
island) the functional dependence of the thermopower on the gate voltage is qualitatively 
changed while its amplitude remains on the same order of magnitude. This is in clear con- 
trast to the corresponding results for the conductance where the most pronounced effect is 
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a suppression of the amplitude with increasing gap size 

Thermopower of single- electron devices - In the following we investigate the regime of 
single-electron tunneling through a superconducting island with a charging energy E c that 
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is large compared to the temperature T, i.e. E c = e 2 / (2C) ^> T, where — e is the electron 
charge and C the capacitance of the island. The temperature is assumed to be larger 
than the crossover temperature for parity effects Q|. The electronic spectrum of the grain 
in the normal-conducting state is assumed to be continuous and the conductances of the 
tunnel barriers are much smaller than e 2 /h. This implies that cotunneling processes can be 
neglected and sequential tunneling dominates. Taking into account the external electrostatic 
potential 4> oc V g imposed by the gate voltage V g , the total electrostatic energy of the island 
can be expressed as 



E n {(j>) = E c (n 2 -2n 



C 



(1) 

where n is the number of excess electrons on the island. To add one electron to the island 
an energy u n {<f>) = E n+ i((f>) — E n (<f>) is required. At low temperatures T <C E c the electronic 



transport is Coulomb-blocked. A current flows only at potentials 
by the condition u n (</>„) = 0. With Eq. ((TJ) one finds 



') n where n is given 



(2) 



The current through the device can be either due to a transport voltage V or a temper- 
ature difference AT = T\ — T r between the two leads (see Fig. 1). Throughout this work we 
consider the linear response regime, i.e., eV/E c <C 1 and AT /T <C 1. 
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FIG. 1: The NSN SET consists of a superconducting island (S) which is coupled to two normal- 
conducting leads (N) via tunnel barriers. The electrostatic potential of the island can be controlled 
by the gate voltage V g . The current through the system is due to the bias voltage V or a temperature 
difference AT = TJ — T r between the two leads. To measure the thermopower S = —V/AT as a 
function of the gate voltage V g the bias V is adjusted in such a way that the corresponding current 
exactly cancels the current due to the temperature difference. 
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For zero temperature difference the linear response to the voltage V is given by the 
conductance Gy. As a function of the potential (f) it shows the well-known CB peaks at 
= (f) n . In the general case with AT 7^ the current is I = GyV + G^AT. The 
thermopower of the system is defined by the ratio of voltage and temperature difference for 
vanishing current, i.e., 



S = — lim 



= p-. (3) 
1=0 u v 



at-*o AT 

Following Matveev [15] a rather intuitive interpretation of the CB oscillations of the 
thermopower can be given in terms of the average energy (£) of the electrons that carry the 
current through the system 

S=-& . (4) 
eT v ; 

Here, we briefly review the argument in the NNN case (normal-conducting island) consid- 
ering a two-state approximation (valid for T <C E c ). If in equilibrium there are n electrons 
on the grain two different transport cycles are possible: one can first add an electron and 
then remove it again [n — > (n + 1) — > n] or vice versa [n — *■ (n — 1) — ► n]. In the first 
case the average energy is given by the difference of charging energies, (£) = (E n+ i — E n )/2, 
while in the second case it is (£) = (E n — E n _ 1 )/2. The potential <fi determines which of 
the processes is more likely to occur. For example, at <p > <p n the probability of having n or 
n + 1 electrons in the grain is higher than the probability for n — 1 electrons. Thus, the first 
process dominates, and according to Eq. (J3J) one finds the well-known sawtooth behavior 

Snnn(0) = for C\<f>- 0„|/e < i , (5) 

see inset Fig. |21 The extrema of the sawtooth at e\<p — <p n \ ~ E c are rounded due to the 
finite temperature. 

Thermopower of NSN SET - In the following we consider the case of a superconducting 
grain with a gap A < E c . In a stationary state the currents through the left and the right 
tunnel junction are equal, i.e., 1 = 1\ = l r . In order to evaluate the thermopower we make 
use of relation Q and calculate the linear response of the current to a voltage V or a 
temperature difference AT = T\ — T r . According to the "orthodox theory" |l( the current 
through the system can be written as 

i = -eJ2Pn [rr^ 1 - rr n+1 ] (6) 



4 



where P n is the stationary probability for finding n electrons on the island, r™^ n_1 is the 
tunneling rate of an electron from the island to the right lead, and T^ n+1 denotes the 
tunneling rate from the right lead to the island. The rates T r and correspondingly Ti can 
be written in terms of the Fermi-function f(x) = 1/[1 + exp(a;)] as 
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(7) 
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(8) 



Here, B is the energy of the quasiparticles in the superconductor, and G r j is the conductance 
of the right and left tunnel junction, respectively (see, e.g., Ref. The first term in 

Eq. (J7J) corresponds to the annihilation of a quasiparticle while the second term yields the 
contribution due to the creation of a quasiparticle. 

The probabilities P n in Eq. (jUJ) can be obtained from the stationary solution of a kinetic 
equation. They obey the relation P 



i n— >n+l 



n+l 



pn+l-+n _|_ pn+1- 



(9) 



with J2 n Pn — 1- The set of Eqs. allows us to calculate the transport coefficients Gy 

and Gt- Together with relation (j^J) the thermopower S can be obtained as a function of the 
potential for different parameters A and T. The results of our calculation are shown in 
Fig. El Compared to the sawtooth-like behavior of the thermopower in the NNN-case (see 
inset Fig. |2J) the shape is significantly changed for gaps A smaller than the charging energy 
E c . For increasing A — > E c the extrema of S(4>) move clearly away from <p n ± e/(2C). 

We emphasize that the qualitative change of the thermopower as a function of the po- 
tential (j) due to a finite gap is much more pronounced than the changes in the conductance 
where one merely finds a broadening of the Coulomb-blockade peaks together with an overall 
exponential suppression of the current. Note however that the order of magnitude of the 
thermopower remains the same independently of the gap size. This behavior can be under- 




FIG. 2: Thermopower S((fr) of the NSN setup with A < E c for different values of the supercon- 
ducting gap A at a temperature T/E c = 0.1. Note that the conductance peak would be located at 
(j) = (fro = 0.5e/C [cf. Eq. l[2*|)]. The inset shows the corresponding result for the NNN SET. 

stood by analyzing Eqs. in the low-temperature regime where the major contribution 

to the electronic transport is due to only two charge states. 

Two-state approximation - For low temperatures T <C E c there are at most two proba- 
bilities that assume a finite value for a given <fr, e.g., Pq and Pi for C\(fr — (fro\/e < 1/2. All 
remaining probabilities are exponentially small in the parameter E c /T. To first order in 
the perturbations V and AT the current © can be expressed in terms of the unperturbed 
probabilities P^ and the exact rates T^ r} Eqs. (J7J) and (jHJ), as 

p (0) (g^ 1 - cvrr 1 ) + (io) 

P/ 0) {GrT™ - G,!*- ) 



The zeroth-order probabilities can be obtained from Eq. and the condition Pq+P^ = 1 
which yields Pq° ~ /[=p/o(0)/T] if exponentially small corrections are neglected. 

First, we discuss the thermopower for potentials close to (fr , i.e., C|0 — o |/ e < A/ (2P C ) 
where |wo(0)| < A- In this range of 0, the changes in the thermopower due to the super- 
conducting properties of the island are most prominent, see Fig. El For gap energies that 
are not too small A < E c we can neglect contributions that are exponentially small in A/T. 
Thus we find from Eqs. (J7J), (JHJ), (fTUj) together with relation Q the asymptotic result 



~ ^r- 1 7-t tanh 

eT \ tt o (0) 



•up(0) 
2T 



(11) 
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Here, A(T) = —A K[(A/T)/Ki(A/T) where K[(x) is the derivative of the Bessel function 
Ki(x). In contrast to the NNN case we find S = not only for iio(0) = but also at 
|<u o (0)| w A, see Fig. El Between these two zeros the thermopower reaches its extrema at 
|w max | ~ 2T _1 arccosh-y/ A/ (2T). This novel behavior that the slope of S((j)) changes its 
sign at = O occurs even for small gap values A ~ IT . In the low-temperature limit 
A/T ^> 1 the temperature-dependent "effective gap" A(T) in Eq. (|TT|) is simply replaced by 
the constant gap A. On the other hand, Eq. (fTT|) also reproduces the limit A — ► correctly 
as it gives A — > T leading to the NNN result of Eq. 

Next we consider the thermopower in the range A < |«o(0)| < E c , i.e., A/(2E C ) < 
C|0 — 0o | /e < 1/2. In this case Eq. (fTUj) and Eq. (JH1) approximately yield 
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(12) 



From this result we find S = for |i*o(0)| = A which agrees with Eq. (|llj) in the low- 
temperature limit where A — >• A. The second term in Eq. (fT2*|) depends only weakly on 
but gives an overall shift by A/2, see Fig. |21 Therefore, the dependence of S on the 
potential is almost linear so that the extrema for Wo(0) ~~ > ±E C are approximately given 
by S max ~ t(E c — A)/(2eT), respectively. Similarly to the NNN case there is a thermal 
smoothing close to the edge of the interval where E c — \uq(4>)\ ^ T. This is because 
for these values of charge states with n = -1 or n = 2 electrons on the island become 
important. For small gaps A — > also Eq. (fT2*)l reproduces the NNN result of Eq. 

Interpretation of results in terms of average energy - The asymptotic results (jllj) and (fT2*|) 
can be intuitively understood in terms of Eq. (@J) by considering the average electron energy 
(£} of the dominating transport mechanism. In Fig. El we schematically present the transport 
mechanisms for potentials such that |«o(0)| < A. The two processes corresponding to the 
rate rj - * , Eq. (J7J), are shown. The Fermi level of the lead is chosen to be zero. The average 
energy of the electrons involved in process A is then (£4) w ^o(0) + A > 0. On the other 
hand, process B involves the breaking of a Cooper pair: one electron tunnels to the lead 
while the other remains as a quasiparticle on the island. Hence, the average energy of the 
outgoing electrons is given by (£b) ~ «o(0) — A < 0. 

As can be seen from Fig. |SJ the two processes are not equally likely to occur due to the 
different occupation numbers in the island and the lead. We can read off a low-temperature 
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FIG. 3: Tunneling processes from the superconducting island to the right lead for |uo(</>)| < A: 
annihilation of a quasiparticle (A) and creation of a quasiparticle (B). The corresponding average 
energies (measured from the Fermi energy of the lead) are (£4) > and < 0. Depending on 
which of the processes is dominant the sign of S results according to Eq. (|H>. 

estimate for the probabilities from Fig. El Pa oc exp[— A/T] and pn oc exp[— (A — u )/T]. 
Thus, relation (jlj together with the condition p& +Pb = 1 leads to the estimate 



S 



-^(Pa{U)+PB{i B )) 



1 

eT 



uq(4>) — A tanh 



up{4>) 

2T 



(13) 



This corresponds precisely to the result (jllj) in the low-temperature limit A/T ^> 1 where 
A(T) -> A. 

If the potential increases further such that Uo((fi) < —A, process A clearly dominates 
as pa ^ Pb an d one can neglect process B entirely. However, the tunneling of low-lying 
quasiparticle excitations is also strongly suppressed. By including the energy dependence of 
the density of states in the superconductor we find an estimate for the average energy (£4) 

(£(!)> 



(U) ~ «O(0) + 



-u (<t>) 



( E (k)) 



dE E 



with 



E 



y/E 2 - A 2 



(14) 
(15) 



Including this estimate into relation (jlj) directly yields the result (|12l) . 

In conclusion, we have developed a theory for the thermopower S of a NSN SET with 
A < E c . The numerically exact results are presented in Fig. El and the asymptotic low- 
temperature behavior is given in Eqs. (jllj) . ()12p . We showed that these results can be 
understood on the basis of the average-energy interpretation Eq. (J3J). We mention that the 
current- volt age characteristics of a NSN SET in the relevant range of the parameters dis- 
cussed in this work have already been studied in experiments [la [l7|. It should be well 
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within reach of present-day nanotechnology to experimentally detect the sensitive depen- 
dence of the thermopower on the gap size summarized in Fig. [21 
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